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^ ! 1. Introduction 

< 

The study of charm electroproduction has become an important issue in the ex- 
traction of parton densities in the proton. The reason is that the charm content 
J> i F2,c{x , Q"^ , m"^) , where m = mcharm, has grown from around one percent of the con- 

^ I tribution to the structure function F2{x,Q'^) in the x and region of the EMC 

CN I experimental to around twenty-five percent in the x and region of the Hli and 

^ ■ ZEUsi experiments at HERA. Therefore the analysis of parton densities in the pro- 

■ ton can no longer be done without considering c-quark production. 



Let us begin with a brief review of some of the technical points in the calculations 
of c-quark electroproduction in QCD. We consider F2^c{x,Q'^,m'^) which is usually 
, much larger than Fl c{x, Q^, m^). We work in a three flavour number scheme (TENS) 

Qj I where the u, d and s are light mass quarks (n/ = 3). The c c pair is produced from 

the gluon by the Bethe-Heitler process (photon-gluon fusion) in leading order (LO). 
In next-to-leading order (NLO) both the Bethe-Heitler and the Compton processes 
^ ■ involve all light mass quarks (and antiquarks) and the gluon. A NLO calculation 

■ organizes the contributions from the various Feynman diagrams according to the 

following formula 



FU3,.,<^^m^)=. r"-(ii:4k3.-.A.=)4,(3.--.^4) (1) 



9 ' 9 



+G(3, -,/i^)L|,(3,., ^, ^) + A(3, -,/i^)Lg(3 

I - S(3, )^2.(3, -, -) + G(3, -, ^^ )if,,(3, - 



where Cc is the charge of the c-quark. The variable z is the partonic longitudinal 
momentum fraction and Zmax = Q"^ /{Q"^ + 4m^). The function A is the non-singlet 
(with respect to the flavour SU(3) group) combination of light-mass parton densities, 
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at the scale /i. The function S is the singlet combination of these densities while 
G is the density for the gluon. Further L2,fc and H2^k{k = q,g) are the singlet and 
non-singlet heavy-quark coefficient functions at scale fi. Since they contain factors of 
as we have explicitly indicated that they depend on n/ = 3. The L2,fc describe the 
reactions where the virtual photon couples to the light quarks (u, d, and s) and the 
i72,fc the reactions where it couples to the c-quark. 

The NLO contributions to Eq.(l) were originally calculated inB yielding single- 
particle inclusive transverse momentum and rapidity distributions of the c-quarki. 
The functions L2^k and H2^k were only available in the form of two-dimensional in- 
tegrals. To speed up the computation of Eq.(l) we made two-dimensional grids of 
values for L2,fc and H2^k together with an interpolation routine ini. The NLO cal- 
culation was repeated in a completely exclusive fashion ini so that one can plot all 
the distributions containing the c-quark, the c antiquark and the additional parton. 
This program incorporates a fragmentation function so that one can calculate distri- 
butions for D mesons. These programs have been used by many authorgsip 
to discuss the sensitivity of the NLO results to changes in parton densities, renor- 
malization/factorization scales, and the mass of the c-quark. Finally the asymptotic 
formulae for L2^k and if2,fc in the limit ^ were calculated The latter 

formulae enable one to discuss the transition from the calculated c-quark rate in 
NLO perturbation theory in the TFNS with nf = 3 to a four flavour number scheme 
(FFNS) with rif = A and a massless c-quarkEJ,!!^. 

Now we outline the scheme used for these NLO calculations. The renormalization 
of the virtual graphs with loops containing light quarks was done in the MS-scheme 
while loops with heavy quarks were subtracted at zero external momentum. This 
scheme was originally proposed inS so that the light quark sector remains the same 
as in a purely massless pQCD calculation. Hence there is a matching condition at 
the scale jj, = m for both Og and for the parton densities. To explain this consider 
QCD with only massless quarks. Then the number of flavours enters via the factors 
n/ in the /3-function and in the gluon splitting function Pgg. We solve the equation 

= (2) 



din n 
where 

/5o = ^(yC^A-^T;n^) , Pi = ^{fc^ - iC^T^nj - fc^Tju^) , (3) 

by introducing the Aj^(n/). The solution of Eq.(2) is 



_ 1 Pi ln(ln(/iVAy ) 

l2 A n3 lv.2/,,2/A2 



+ - (4) 



A matching condition is therefore necessary to go from a TFNS with rij = 3 to a FFNS 
with Hf = 4. We have to change the value of as we cross the scale = m to keep 



the condition that a<((4, = m) = 0^(3, /i = m). However heavy particles cannot be 
ignored in QCD because they exist in vacuum fluctuations. It is well-knownS that 
effects from very heavy particles can be absorbed into an equivalent field theory at a 
smaller scale by an appropriate renormalization of the coupling constants and fields 
(decoupling of heavy flavours). QCD is an unbroken gauge theory which does not 
have heavy particle decoupling when renormalized in the MS-scheme as all quarks 
have m = 0. However the top quark cannot be considered massless and virtual top- 
quark loops yield explicit terms in In/x^/m^ (here m = mtop) so there has to be a 
matching condition between say a five flavour theory and a six flavour theory. This 
two-loop matching condition on ag was first worked out inli and corrected in@. It 
reads 

Q«(6,/i) ^ as(5,/i) ^ / a,(5,/x) N2l^ \n(—] + ( ^^^^^ll^Y \-T^ (1^] 
+1(5C^T, - 3C,T,) In (^) + f^T.C, - ^T,C^] + ... . (5) 



Note that as(6, = m) = as{5, = m) in order but not in order a^. 

The parton densities must satisfy similar matching conditions to Eq. (5). In a 
NLO calculation the collinear singularities are regularized by n-dimensional regular- 
ization. We set the mass factorization scale equal to the renormalization scale /i for 
simplicity. Only the gluon-gluon splitting function 

Pgg{nf, z) = 2Ca[z{-^^^ + + z{l - z)] + -{lie A - ^Tjnf)5{l - z) , (6) 

depends on the number of light quarks Uf. The c-quark is considered heavy with 
mass m = mcharm and we calculate the L2^k and H2^k in Eq- (1) in the TENS. There 
is no charm density at the scale /i = m but it does exist for /z > m. This means that 
the light flavour densities /^(a;, /i), k = u,d, s, it, d, s satisfy 

/fc(4,a;,/i = m) = /fc(3,a;,/i = m) . (7) 

One can examine the DGLAP equationslll for the parton densities when one chooses 
a scale with fi > m. By taking the difference of the solutions with Uf = 4 and Uf = 3 
and using Eq. (6) one finds that only the gluon density changes in order a^, namely 



1 

/g(4, X, /i^) = fg{3, X, vr?) 1 - — r/as(3, /i) In — , 



and there is a c-quark density, proportional to q;s(3, /i) ln(/i^/m^), which grows at the 
expense of the gluon density. The changes in the light flavour densities begin in order 
a^. Hence at scales ^ m we can switch to a EENS which contains a (massless) 
c-quark density and a different momentum sum rule. Eq.(8) is a matching condition 
on the gluon density between the TENS and the EENS at the scale fi = m. 
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Figure 1: Scale dependence of F2^c{x = 0.05, Q = lOGeV/c). 

A way to implement a smooth matching of F2^c{x, Q^, m^) from the TFNS in Eq. 
(1) to a charm density description in the FFNS was proposed ini and we call it the 
ACOT scheme. Their scale choice, which we will use in the rest of this article, is 

= + kQ\l-myQ'')'' for > m\ 

= for < (9) 

with k = 0.5, n = 2 and m = 1.5 (GeV/c^). They proposed that F^^^^'^^\x,Q'^ ,m^) 
should consist of three terms. The first one is the TFNS LO process in Eq.(l), which 
has the correct kinematics near threshold. It contains the convolution of the gluon 
density in the proton with the LO hard scattering cross section. The third term is 
a FFNS charm density, i.e., the coefficient function in the proton structure function 
F2{x, Q"^) is simply a (5-function, and this should be the best description at ji^ ^ rn?. 
Finally the second term is the product of two convolutions. Namely that of the 
gluon density in the proton with the massless c-density in the gluon and with the 
LO hard scattering term (a 5-function). This second term is taken with a negative 
sign so that it cancels the FFNS c-density term at small scales /x^ ~ m? while it 
tends to cancel the TFNS term when /i^ ^ m? . The sum of these three terms is 
more stable under a variation of the scale /i than each term separately. In Fig.l 
we show the scale dependence of F2^c{x-,Q'^ iTn?) for the TFNS LO term from Eq.(l) 
(labelled EXACT,(1)), the ACOT description (labelled AC0T,(1) and the FFNS 
charm density (labelled PDF,(1)). We choose the same CTEQ parton densities^ as 
in@. The monotonic decrease in the EXACT result, which is due to the decrease in 
is modified in the other two approaches. In Fig. 2 we plot the three results for 
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Figure 2: F2 c(a: = 0.01) as a function of Q. 



F2,c{x = 0.01, Q^) versus Q = \/Q^. One sees that there is a large difference between 
the EXACT result (TFNS) and the PDF charm density result (FFNS) especially as 
Q increases. However the charm density result cannot be used at small Q where it 
is negative. As expected the ACOT result interpolates between the EXACT result 
at small Q and the charm density PDF result at large Q. The ACOT formula has 
been incorporated into recent parton densities (the set CTEQ4HQli). Also there 
are now three light flavour (CTEQ4F3) and four light flavour (CTEQ4F4) global fits 
to parton densities. The former is basically equivalent to the parton density set in 
GRV94i. An examination of the CTEQ4HQ set shows that the gluon is diminished 
from previous setJ^ll and that all the light quark densities have changed, not just the 
one for charm. This is due to the interplay between the gluon and the light flavour 
densities in making a global fit to all the experimental data. 

As the ACOT description contains algebraic terms in m, it is not obvious how 
to generalize it to higher orders in a^. (One proposal has been made in @ which we 
cannot discuss here for lack of space). One would like to retain as much as possible a 
pQCD description based on massless (3 and 7 functions in the renormalization group 
equation, and work in the MS scheme. Inlll we have worked out the two-loop matching 
conditions on the flavour densities by the following procedure. Analytic formulae for 
the functions if2,fc in Eq-(l) are not known. In the limit that ^ they reduce to 
expressions like a{z) In^ (Q'^ /m^)+b{z) ln((5^/m^)+c(z) (where we call them if^f^'^^). 
If we can find a{z), b{z) and c{z) then we can mass factorize the logarithms into terms 
involving ln(Q^//i^_) and ln(/x^/m^). If we then add the result to the two-loop light- 
flavour corrections^^ to the proton structure function ^2(8, x, Q^, /i^) then the terms 
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Figure 3: Scale dependence of F2^cix — 0.05, Q — lOGeV/c). 

in ln((5^//U^) should rccombine to yield the two-loop corrections to F2{4:,x,Q^, fi'^). 
Note that F2,c niTist be totally inclusive so we have to add the light quark vertex 
corrections where the gluon propagator contains a charm quark loop, which yields 
terms in S{1 — z). When the correct pieces have been identified we are left over 
with terms in ln^(/x^/m^) and ln(/x^/m^) which cannot be absorbed into the known 
matching conditions on as- We can then define two-loop matching conditions for the 
parton densities to incorporate these pieces, analogous to the two loop relations for 
as in Eq. (5). We actually derived L^l"^^^, i/^f^^'^ and the relations between the 
parton densities by calculating all two-loop matrix elements of the operators in the 
OPE, where one loop contains a massless quark and the other loop a massive quark. 

The result is that above the scale /i — m we have new FENS parton densities 
{uf — A) which are calculable from a TENS set {uf — 3). All densities are changed. 
The light quark densities pick up additional terms in order a^, which are numerically 
rather small, but are required for the EENS parton densities to satisfy the momentum 
sum rule. The EENS charm density is given in terms of the previous TENS densities 
and involves terms in order as and a^, i.e., 

/J™S(4, X, = /4(4, X, f,') + /4(4, X, f,') 
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Figure 4: F2^c{x = 0.01) as functions of Q. 
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where 



4PS,(2) 
^cq 



+ 



m 



[1 + z) \n^z- 



16 ^ ^ 16 2 

Ul + z)\nz 4 + 4z H z^ 

^ ' 3z 3 



In^ 



64 



160 



m 



448 



40;z + —z"^ In 2 - — + 16 - 48;z + 

3 9z 9 



32Si,2(l - ;z) + 16 ln;zLi2(l - z) - 16C(2) Inz - - In^ ^ 

3 



+ 



2 

'32 



32 o\ . s / 32 32 n\ , / 16 o\ , 2 

-8z - —z" \ Li2(l - z) -f I -— - 8 + 8z + —z^ j C(2) + I 2 + lOz + —z^ j In^ z 



56 88 448 

\ z H 

3 3 9 



z \\n.z 



448 



124 



1600 



-z + 

21z 3 3 27 



(11) 



The other contributions to Eq. (10) and the relations between the other parton 
densities are available inS Eq. (11) shows that the charm density has a discontinuity 
at the scale n = m when one considers order terms, just as in Eq. (5). 

The above formula is for a FFNS renormalized in the MS prescription. Using the 
FFNS one can now calculate -^^^^'''^''(4, x, Q^, m^) in the region where ^ and 
include the higher order QCDc^. We expect this description to be valid in the region 

^ where the TENS formula in Eq. (1) is dominated by the mass factorization 
logarithms. We have checked numerically that the asymptotic result for the NLO 
formulae, which we call F2^^^^''^'^\3, x, Q^, m?) is a good approximation to the exact 
NLO result in Eq.(l), which we call F^^^'""'"'^^''(3, a;, Q^, m?) when ^ 10 m^. Hence 



rr[^= 1 .5 GeV/ c 



x = 0.0005(X4000) 



10- 



10' 



10 



1 



10 



-2 




Q^(GeV=) 



1 



Figure 5: The functions F2^^'^'^\4:, x, Q^) dashed hne, F^™^'*'^''(a;, Q^, m?) sohd hne 

and ^^^^^"""'^^^(3, x, Q^, m?) dotted hne, plotted as a functions of for fixed x. The 
data points are from the EMC, HI and ZEUS experiments. 



above ~ 10 we have the option of retaining the TFNS description for charm 
or to switching to the FFNS which sums large logarithms in ln(Q^/ /i^) . 

It now remains to find a formula which allows us to interpolate between the TFNS 
and the FFNS and is valid in all orders in a,. Ini and0 we have proposed the variable 
flavour number scheme (VFNS) formula 

Q^ m^) = Fj,^^^^(3, x, Q\ m^) -^^^^^^(3, x, Q^ m^) + F™^(4, x, Q') .(12) 

This formula is designed like the ACOT one. We have all the pieces to test it in 
order a^. (PuU details are available in0). The TFNS ASYMP term cancels the TFNS 
EXACT term at large leaving the FFNS charm density PDF term. However neither 
the ASYMP term nor the charm density PDF term have the correct kinematics (they 
also have different Uf and electric charges) to cancel each other exactly at small (near 
threshold). Using as input the same CTEQ densities^ as in§, we have calculated 
(1) i7|^^ACT,(2)^ ^p^g formulae in Eq.(l) (2) ^2^™^'^^^ the VFNS formula in 

PDP ('2) 

Eq.(12) and (3) ^ , the massless FFNS PDF result. In (3) we have evaluated 
the two convolutions to go from a TFNS to a FFNS and then to add the higher order 
correct ionsii. In Figs. 3 and 4 we show the analogous results to those in Figs. 1 
and 2. One sees that the scale sensitivity is mostly reduced for the TFNS EXACT 
result. The PDF result is worse at small Q, which is a consequence of the fact that 
some of the higher order contributions in QCD are negative. The VFNS result Eq. 
(12) interpolates nicely between the other two. Finally in Fig. 5 we compare the 
three descriptions for Fo^ with the data from the EMC, HI and ZEUS experiments. 
Here we use the GRV92B3 parton densities. One sees that the VFNS interpolation is 
quite smooth as compared to that i n@. The massless FFNS results fall dramatically 
at small Q, due to the choice of /i in Eq. (9). All three curves are within the present 
experimental errors of the HI and ZEUS experiments. Therefore we will have to wait 
for more precise data before we can make a firm statement about which theoretical 
approach is the best. 
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